Elasticity of Stiff Biopolymers 
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We present a statistical mechanical study of stiff polymers, motivated by experiments on actin fil- 
aments and the considerable current interest in polymer networks. We obtain simple, approximate 
analytical forms for the force-extension relations and compare these with numerical treatments. 
We note the important role of boundary conditions in determining force-extension relations. The 
theoretical predictions presented here can be tested against single molecule experiments on neu- 
rofilaments and cytoskeletal filaments like actin and microtubules. Our work is motivated by the 
buckling of the cytoskeleton of a cell under compression, a phenomenon of interest to biology. 
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In recent years, statistical mechanics of semiflexible 
polymers has emerged as an active area of research. This 
has been triggered by single molecule experiments de- 
signed to understand the role of elasticity of these poly- 
mers. Elastic properties of polymers are of importance in 
biology as in the structure of the cytoskeleton a biopoly- 
mer network which controls cell mechanics [J Q . The pa- 
rameter which determines the stiffness of a polymer is /9, 
the ratio of its contour length L to the persistence length 
Lp. While the entire range of (3 is of biological interest, 
in this paper we focus attention on rigid filaments such 
as actin filaments and microtubules which constitute the 
cytoskeletal structure and serve as tracks for motor pro- 
teins like myosin and kinesin[3, Recently, filaments of 
intermediate rigidity like neurofilaments have also been 
studied in some detail It has been shown that some 
remarkable features of single stiff filament bending re- 
sponse are relevant to crosslinked biofilament networks 
[21 ■ A good understanding of the elastic properties of 
biopolymers at the single molecule level, is essential to a 
study of polymer networks. 

There are two classes of experiments which probe 
the elasticity of single biopolymers. In one class of 
experiments [5[ a semiflexible polymer molecule is pulled 
and stretched to study its "equation of state" by mea- 
suring its extension as a function of applied force. In the 
other class of experiments, one tags the ends with fluores- 
cent dye[^[3| to determine the distribution of end-to-end 
distances. Such experimental studies provide valuable in- 
sight into the mechanical properties of semiflexible poly- 
mers. A good theoretical model is needed to correctly in- 
terpret these experiments. A simple and popular model 
which captures much of the physics is the Worm Like 
Chain modelfi]. 

In this paper we analyze the bending degrees of free- 
dom of a stiff polymer where at least one end of the 
polymer is clamped. By this it is meant that the tangent 
vector at this end is kept in a fixed direction. Just as in 
the case of a stretched polymer[9i|, the tangent vector of a 
stiff polymer executes small wanderings around this fixed 
direction. The theoretical analysis for the statistical me- 



chanics of a stiff polymer, clamped at least at one end 
is similar to that of a polymer in the high stretch limit. 
We refer to this approximation as the paraxial approx- 
imation. This approximation has been previously used 
to study the elasticity of twist storing flexible stretched 
polymers [I^, [111, Hi in the paraxial worm like chain 
model (PWLC model). In an earlier paper ^9j| we had 
an "exact" numerical scheme for a semiflexible polymer 
with free boundary conditions for the end tangent vec- 
tors. The main difference in the numerical scheme pre- 
sented here is that we impose boundary conditions on the 
tangent vectors at the ends. As a general rule, in the stiff 
regime the convergence of the numerical scheme is poorer 
and therefore one needs to use larger matrix sizes. The 
search for simple analytic forms to describe the elastic 
properties of stiff polymers is therefore well motivated. 

For stiff polymers, the experimentally measured mean 
values crucially depend [3| on the precise choice of the 
ensemble. This is due to finite size fluctuation effects, 
which are entirely absent in the elasticity of a classical 
rod. For instance, one gets qualitatively distinct features 
in force-extension curves depending on whether the force 
or the extension is held constant in an experimental setup 
[isl . [igI . 17 1 . This is an aspect of stiff polymer statis- 
tical mechanics which is both theoretically challenging 
and experimentally significant. In this paper, we remain 
throughout in the Gibbs ensemble, where the applied 
force is held fixed and we measure the mean extension. 

The organization of this paper is as follows. We first 
present results based on an "exact" numerical scheme for 
stiff polymers for two different boundary conditions, one 
in which both ends are clamped and the other in which 
one end is clamped and the other end is free. We then 
present simple analytical forms for these two cases. Our 
main results are displayed in Figs. 1-3 comparing the 
numerical scheme with the analytic formulae for force- 
extension relations. Finally, we end the paper with a 
concluding discussion of the buckling of stiff polymers 
and the consequent breakdown of the paraxial approxi- 
mation. 

Our starting point is the Worm Like Chain (WLC) 
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model in which the configuration C of the polymer is de- 
scribed by a space curve x{s), with s the arc-length pa- 
rameter (0 < s < i) ranging from to L, the contour 
length of the polymer. The tangent vector t ~ dx/ds to 
the curve is a unit vector 

i.i = 1 (1) 

and the curvature of the polymer is given by k = \dt/ds\. 

One can study the case of stiff polymers using a combi- 
nation of analytical and numerical techniques [9| . Let one 
end of the polymer be fixed at the origin and a stretching 
force F be applied in the z direction, which we refer to 
as the north pole of the sphere of directions. Introducing 
a dimensionless force variable / — where fc^T is 

the thermal energy we can express the partition function 
Zif) as 

Z{f) =M j V[i{T)]e' r dr[l/2(di/drf-fh] (3) 

where (3 = L/ Lp. Eq. ^ can be interpreted as the path 
integral representation for the kernel of a quantum par- 
ticle on the surface of a sphere at inverse temperature 
(3. Thus we can express Z[f) as the quantum amplitude 
to go from an initial tangent vector Ia to a final tan- 
gent vector Ib in imaginary time (3 in the presence of an 
external potential —fcosd: 

Z{f, tA, ie) iA\ exp[-/3F/]|fB > (3) 

The Hamiltonian Hf — /cos 6* is that of a rigid 

rotor 18] in a potential. In the absence of a force, the 
free Hamiltonian is Hq — — ^V^. By choosing a stan- 
dard basis in which Hq is diagonal we find that H is 
a symmetric tridiagonal matrix with diagonal elements 
Hi,i = l{l + l)/2 and superdiagonal elements — 
f{l + 1)a/1/((2/ + l){2l + 3)). Inserting a complete set 
of eigenstates of the free Hamiltonian into Eq. ([3]), we 
find Z{f,iA,iB) = 

< tAlipn >< 'fpnl exp[-f3Hf]\i{jm >< V-mlis > (4) 

Tn,n 

where M-^ = exp[—(3Hf]. From this general form, we 
can compute the partition function in the present cases 
of interest. 

In Rcf . ^ we had studied the elastic properties of poly- 
mers with free boundary conditions: the directions of the 
tangent vectors at both ends were integrated over. In the 
present paper, we will fix the tangent vector at the ends 
(one or both) to lie along the z direction. To implement 
this numerically, we have to evaluate the eigenfunctions 
in Eq. ([4]) at this value of t. 

{i)both ends clamped: Ia = is = z. While a complete 
set of eigenstates are labelled by {l,m), only the m = 



terms contribute here because of azimuthal symmetry 
and we have 

Z{f, z,z)^Yl UiM[,, Uv = U.M.U (5) 
1,1' 

where Ui = ^[^. 

[a) one end clamped: Integrating Eq. ([4]) over is, we 
find that 

Z{L tA) = f^'^^/o = iU.Mf )^ (6) 
I 

Both Eqs. ^ and ^ are suitable for numerical im- 
plementation. Hf is an infinite symmetric matrix. We 
truncate it to finite order N and choose TV large enough 
to attain the desired accuracy [l9| . 

While this numerical method is effective, it has a lim- 
itation in describing stiff polymers due to the poor con- 
vergence of statistical sums in Eq. Q. For stiff poly- 
mers a convenient and accurate analytical approximation 
scheme can be developed as shown below. 

For a stiff polymer with one end clamped along the z 
direction, we can approximate the sphere of directions 
by a tangent plane at the north pole of the sphere as the 
angular coordinate 9 always remains small. Introducing 
Cartesian coordinates i^i = 9 cos cf) and ^2 ~ 6 sin cj) on 
the tangent plane at the north pole one can express 
the small 9 Hamiltonian H a.s H = Hp — f where Hp is 

^p = ^P6' + ^P6' + |(ei' + 6'); (7) 

Notice that Hp is the Hamiltonian of a two-dimensional 
harmonic oscillator with a frequency uj — \fj . For a sin- 
gle oscillator in real time the propagator is given bv[20|: 

F{T) ^;^[(^' + cos^?^ - 2^.0] (8) 

where F(T) = y=ij=^. 

For the sake of convenience we set Lp = 1 so that 
/3 = L. 

Case (i): both end tangent vectors clamped along the z- 
direction: Setting = = in Eq. ([8]) and continuing 
the expression to imaginary time results in the trigono- 
metric functions being replaced by hyperbolic ones. We 
can express the partition function Z(f) as exp (/?/) times 
the product of the propagators of two independent har- 
monic oscillators: 

^(/) = y7exp(/3/)/(27rsinh(/3y7)). (9) 

in Euclidean time /3; the free energy is G{f) = 
-logZ(/)//3 

= l^^"9f -f + 4^05(277) + ^log[smh{P^)] (10) 
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The mean extension < ( >=< z > /L = —dG( f)/df is 
given by (See Fig. 1). 

< C >= 1 + 1/(2/3/) - coth(/3/7)/(2/7). (11) 

where < C > is the z component of the extension (or 
the end-to-end distance vector). Note in Fig. 1 that the 
analytical form agrees with the numerical scheme to an 
accuracy of about 1%. 

Case (a): A stiff polymer with one end tangent vector 
pointing along the z-direction and the other end free: 

In this case the propagator for the harmonic oscillator 
has to be integrated over the final coordinates ^/ and 
evaluated at = 0. The partition function in this case 
turns out to be 



Z(/) = exp(/3/)/cosh(/3y7). 



(12) 



From the expression of the partition function, we get the 
free energy 

G(/) = -/ + ^/o5[cosh(/3/7)] (13) 

and differentiate it with respect to / to get the force - 
extension relation - 



< C >= 1 - tanh(/3V7)/(2 V7)- 



(14) 



Note that even at zero force, there is a nonzero exten- 
sion, because of the boundary condition and the stiffness 
of the polymer. Fig. 2 shows a comparison of the force 
extension curves for the two boundary conditions. We 
find, as expected, that for the same force, the extension 
is larger for the more constrained boundary condition 
[Case(i)] compared to a less constrained one [Case(ii)]. 
For positive forces the paraxial approximation is very 
good and the forms are displayed in Eqs. (fTTjl and (fT4l) 
For large positive forces they become < C 1 

[l3| as expected. 

For negative forces, the hyperbolic functions appearing 
in Eqs. ^ and (|12p go over to circular functions. For 
instance, for the case in which both end tangent vectors 
are clamped along the z direction, for negative /, our 
simple analytical form for the partition function reads 



1 

2^7 



zif) = 



(15) 



27rsin(/3^/^) 

and varies continuously with / as / ranges from positive 
to negative values. As the compressive force is increased, 
we find that at a critical value of the force / the extension 
< C > spontaneously decreases. This is the analogue here 
of the classical Euler buckling instability which occurs for 
rods. 

Consider the mean extension versus force relation 
[Eq. fTT]) ] for negative forces (i.e. for compressive forces). 
For negative values of forces Eg. pTj) reduces to 



which can be rewritten in the form 



< C >= l + Pu{x). 



(17) 



where x = fi\[-~] and 



u{x) 



cot(a::) 
2x 



1 

2^ 



The criterion for the onset of the buckling instability 
is the divergence of i9 < ^ > From Ea. (fT7|) this 

is equivalent to the divergence of du/dx^ which takes 
place at a value of Xc — tt. This gives us the following 
expression for the critical force for buckling 2l| : 



< C 1 + 1/(2/5/) - cot(/3v/-/)/(2v/-/). (16) 



Because of the quadratic dependence, the compressive 
force needed to buckle a polymer rises sharply with stiff- 
ness. The mean extension versus force curves displayed 
in Fig. 3 demonstrate the phenomenon of buckling. As 
expected, we notice that as /3 goes up, the magnitude of 
the critical force fc needed to buckle the polymer goes 
down. 

A stiff polymer is energy dominated and its buckling is 
very similar to that of a classical rod subject to identical 
boundary conditions and a compressive force. The ef- 
fect of thermal fluctuations is to slightly "round off" the 
transition from the straight to the buckled configuration. 
This is due to thermally activated processes that permit 
the polymer to overcome the elastic energy barrier. As 
a result the buckling force for a stiff polymer is slightly 
smaller in magnitude than the fc given above. 

There is a long history of the use of path inte- 
grals in the study of polymers [l^, [l^. Such methods 
have been used in the study of elasticity of semiflexible 
polymers 0, 113, IH, [2^. This connection between path 
integrals in quantum mechanics and statistical mechan- 
ics of polymers enables us to import ideas back and forth 
between these two distinct domains. The main point of 
this paper is that standard results in path integrals give 
us new results for stiff biopolymers. Our main results are 
contained in the analytic forms displayed in Eqs. (Pl- fH)) 
and Figs. (1-3). 

In this paper we have theoretically studied the elas- 
ticity of stiff biopolymers. We have studied some cases 
with boundary conditions realizable in single molecule 
experiments. By attaching a magnetic bead to an end of 
the polymer, one can apply forces using magnetic field 
gradients and torques using magnetic fields. By such 
techniques one can impose a variety of boundary condi- 
tions on the polymer including the ones discussed here. 
Recent studies have shown [2] that the elastic behavior 
of such a biopolymer at the single molecular level affects 
the elastic properties of a biopolymer network. This is 
much like the way the structural stability of a roof is de- 
termined by the rigidity of its rafters. In a cytoskeletal 
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structure the end tangent vectors of the stiff biopolymers 
that make up the structure are pinned. A cytoskeleton 
can be viewed as a rephca of a large number N of semi- 
flexible polymers. By studying the elastic properties of a 
single polymer constituting such a network, we can draw 
conclusions regarding the stability of the N polymer cy- 
toskeletal structure. Here we have presented closed form 
simple analytical expressions for force-extension relations 
for a single stiff filament which can be tested against sin- 
gle molecule experiments. These analytical results are 
new and are expected to shed light on the structural sta- 
bility of the N polymer cytoskeletal structure. 

We have also considered the case in which one end 
of a stiff polymer is clamped and the other end is free. 
This is a boundary condition that is more natural to an 
experimental setup for measuring the end-to-end distance 
distribution P(C) of a polymer via imaging of a polymer 
tagged with fluorescent dye. In fact one can construct a 
force-extension curve from the experimental data of P{C) 
versus (. We have theoretically analyzed this case and 
made predictions for experiments in this case as well. As 
in the earlier case, in this case also we have a simple 
analytic form. This is another new result. 

In future, we would like to investigate buckling of stiff 
filaments like actin in greater detail . This is an issue that 
is of relevance at the single molecular level as well as at 
the level of a biopolymer network like the cytoskeletal 
structure and is expected to shed light on its structural 
stability^, and collapse under stress. The stiffness 
and collapse of the cytoskeletal structure of a red blood 
cell "ST] has a direct connection to its functional aspects 
and used for instance, as a diagnostic for detection of 
sickle cell anaemia. In studying the cytoskeletal structure 
it would be most useful to have a good understanding 
of the individual polymers that make up the structure. 
Simple analytic forms give valuable insight into a prob- 
lem and we expect the analytic results presented here to 
provide some fresh impetus to this rapidly growing field 
of semiflexible polymers. 
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FIG. 2: The figure demonstrates that force-extension rela- 
tions depend on the boundary conditions. Note that as ex- 
pected, for a given force, the extension is greater for the case 
where both ends are clamped in the z direction. 
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FIG. 3: Figure shows buckling i.e. spontaneous decrease 
in extension under a compressive force for a stifT polymer 
with end tangent vectors clamped for /? = .5 and fi = .6. 
Note that buckling takes place at a smaller magnitude of the 
compressive force / for a larger /?. 



